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We have measured the in-plane angular variation of nucleation and annihilation fields of a multi- 
domain magnetic single dot with a micro-SQUID. The dots are Fe/Mo(f fO) self-assembled in UHV, 
with sub-micron size and a hexagonal shape. The angular variations were quantitatively reproduced 
by micromagnetic simulations. Discontinuities in the variations are observed, and shown to result 
from bifurcations related to the interplay of the non-uniform magnetization state with the shape of 
the dot. 
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Coherent rotation of magnetization is the simplest 
model of magnetization reversal, proposed by Stoner 
and Wohlfarth in 19480. Coherent rotation predicts 
the value of the switching field H svll of a single-domain 
system as a function of the direction of the external 
field .ffext- For a two-dimensional system with uniax- 
ial anisotropy the polar plot H sv ,i(ip) falls on the well- 
known astroidQ. The full experimental proof for coher- 
ent rotation was given only recently, when nanoparticles 
of high quality and of size small enough to roughly satisfy 
the hypothesis of uniform magnetization could be inves- 
tigated individually 0, 0- Starting from this proof, it 
is now a challenge to understand magnetization reversal 
in increasingly large (and thus complex) systems. The 
simplest ingredient to add to coherent rotation is to al- 
low minor deviations from strictly homogeneous magne- 
tization. The consequences on magnctizationprocesses 
were addressed by numerical micromagnctics|S|, investi- 
gated analytically^ and checked experimentallyQ. The 
next step is now to tackle quantitatively more strongly 
non-uniform systems, those that may display magnetic 
domains and domain walls 0. In such systems a switch- 
ing field i/swi is not the signature of the full reversal of 
magnetization, but instead reflects events like nucleation, 
propagation and annihilation!!!. Few and only partial 
experimental^ and numerical|l0j reports are found on 
this issue. A more detailed study would open the door to 
understanding microscopic details of magnetization re- 
versal processes in macroscopic materials. In this Let- 
ter we present such a study in a model system: sub- 
micrometer Fe faceted dots self-assembled in UHV, that 
have a high structural quality and display simple multi- 
domain states^]! 0|. The angular dependence of the 
i/swi's of a single dot was studied with the micro-SQUID 
techniaue|l3l|. This can be seen as the first experimental 
generalization of astroids for an individual multi-domain 
system. A striking feature is the occurrence of disconti- 



nuities (hereafter named jumps) in H SW i(ip) plots. These 
jumps were reproduced and understood with the help of 
numerical micromagnetism. They result from bifurca- 
tions, related to the interplay of the non-uniform magne- 
tization with the shape of the dot. This also shows that 
a complex _ff sw i behavior does not necessarily result from 
defects, whose role may have been overestimated in the 
literature of magnetic dots made by lithography. 

The Fe(llO) epitaxial dots were fabricated with pulsed 
laser deposition in ultra-high vacuum by self-assembly on 
Mo(110)[8nm]/Al 2 3 (1120). The dots display the shape 
of ingots with atomically-flat facets, bulk lattice parame- 
ter and bulk cubic magneto-crystalline anisotropy K\ fa- 
voring < 100 > axes, however of magnitude much smaller 
than h unM?^n\. The inter-dot dipolar fields are negligi- 
ble with respect to -ff sw i. The remanent state consists of 
flux-closure domains, resulting from demagnetizing fields 
within each dot|l2|. Such domains can occur due to size 
of the dots being well above the exchange length A ex . 0, 
llfij . The in-plane i/ S wi's of a single Fe dot were mea- 
sured below 4K using the micro-SQUID tcchnique^j|. 
For these measurements, the dots were covered in UHV 
by Mo[2nm], followed by Al[2nm] (then 12 hours air- 
oxidized), and a Si[2 nm]\Nb[15 nm]\Si[2 nm] tri-layer. 
Arrays of square micro-SQUIDs with edge lxl /im were 
patterned by e-beam lithography and SFg reactive ion 
etching of the tri-layer. The oxidized Al layer prevents 
ferromagnetic-superconductor proximity effects between 
the dots and the micro-SQUIDs. Although the dots are 
randomly distributed on the surface, their large number 
yields a significant probability to find one suitably cou- 
pled to a micro-SQUID. The location and shape of the 
single dots under investigation were checked a posteri- 
ori by AFM. The size of the dot selected here (FIG. QJ,) 
is 420 x 200 x 30 nm (FIG. GJ)). Micromagnetic simu- 
lations were performed for K (no thermal activation) 
using custom-developed codes, either based on integrat- 
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ing the LLG e quat ion in a finite differences code (rect- 
angular prisms) [lfj or on energy minimization in a finite 
elements code (tetrahedra) . The applied field was in- 
creased step-wise in hysteresis loops. In finite differences 
the sample was divided into cells with uniform lateral and 
vertical size A x = A y = 4.70 nm and A z = 3.75 nm, re- 
spectively. For finite elements 83310 tetrahedra of irreg- 
ular but similar shape were used, with a maximum (resp. 
minimum) volume of 42.29 nm 3 (resp. 12.50 nm 3 ). We 
set Ki = 4.8 x 10 4 J.m" 3 , A = 2 x 10" 11 J.m" 1 and 
M a = 1.73 x 10~ 6 A.ixT 1 in the calculation. 

In the following we call tp the angle between the 
in-plane H cxt and the in-plane long axis of the dot 
[001] (FIG.QJ)). Due to a shape effect, in-plane [lTO] (if = 
90°) is a magnetically-harder direction than [001]. The 
insets of FIG. ^ show micro-SQUID hysteresis loops for 
two angles (tp = +6; +90°). Such loops with negligible 
remanence although with significant hysteresis, are char- 
acteristic of multidomain systems with a limited number 
of domains. Starting from positive saturation the first 
H swl , named hereafter H nuc , is expected to reveal a nucle- 
ation event, e.g. the entry of a magnetic vortexfljj in the 
dot. The second H sw i, occurring at negative fields and 
named H ann , is expected to reveal an annihilation event, 
i.e. the expulsion from the dot of a previously-nucleated 
vortex or wall. FIG.^ shows the experimental angular 
variation H nuc (tp) and H ann (ip). The two-fold symme- 
try results from the elongated shape of the dot. Two 
striking features are observed, that shall be explained in 
the course of the discussion. First, jumps of both H nuc 
and -ffann occur at some angles. Second, depending on 
the range of angles, one or two H nuc and/or H anil are 
observed. 

The jumps of -ff sw i can be understood qualitatively by 
simple arguments. Let us sketch in a quasistatic picture 
the evolution of magnetization M(r) close to an edge 
during the first stages of a hysteresis loop (FIG. EJi-b). 
Starting from saturation, upon decrease of H ext the rel- 
ative importance of the dipolar energy increases. As 
a result M progressively rotates towards the edge to re- 
duce surface charges, and thus reduce The direction 
of rotation, clockwise or anticlockwise, depends on the 
initial direction of M with respect to the normal to the 
edge (imposed by the direction of H cxt ), due to the torque 
exerted by the dipolar field on M. With this picture 
at least two different slightly inhomogeneous magneti- 
zation states, so-called 'leaf states' 0, are expected to 
appear upon decrease of H cxt , e.g. when starting from 
saturation along tp = 0° or tp = 90° (FIG. EM). Bi- 
furcation must occur for at least one intermediate angle 
between these two paths. Then, it is obvious that for 
H C xt applied on either side of this angle, the magnetiza- 
tion pattern will evolve towards very different states, each 
characterized by a different entry point for vortices-and 
thus of edge orientation, explaining a jump in H nuc (tp). 
These ideas were confirmed by micromagnetic simulation. 



In a first attempt the simulations were performed on a 
dot with vertical facets and symmetric ends. FIG. [SJi- 
b shows the static magnetization states just before and 
after H nuc . For ip = 40°, before H nuc the state belongs 
to the class sketched in FIG. Et, as expected. Then two 
regions of strongly non-uniform magnetization develop 
simultaneously, ending up in the entry of two vortices 
at H nuc . As H cxt is further decreased the two vortices 
with opposite circulation move towards the inner part of 
the dot, ending up in a diamond state (FIG. EK)- For 
ip = 50 the state before nucleation looks like that in 
FIG. Eli- The loci of the entering vortices are thus mod- 
ified with respect to the above situation, explaining the 
jump of H nuc , but ending as well in a diamond state, i.e. 
with two vortices (FIG.EtO- 

The above arguments explain the jumps of H nVLC , but 
fail to explain (1) the existence of either one or two H nuc 
and i?ann for some angles (FIG.^) (2) the experimental 
observation of both diamond and Landau states[l2j. i.e. 
with two or one wall or vortex. Indeed in the simulations 
for any tp two vortices appear simultaneously on opposite 
loci, as the dot was assumed to be perfectly symmetric. 
Due to the large dot size these vortices interact weakly 
with each other, thus both enter the dot, ending up in a 
diamond state. In order to refine our interpretations, we 
now report simulations performed on dots with a slightly 
asymmetric shape, similar to that of the AFM observa- 
tion of the measured dot (FIG. QJi). Opposite loci are 
no more equivalent due to the point-reversal symmetry 
breaking. Two situations occur. If opposite loci are sim- 
ilar two vortices still enter the dot, one slightly before 
the other in terms of H exi , however still ending in the 
diamond state (e.g. for tp = +50 °,FIG. EtO- If opposite 
loci arc significantly different, then one of the vortices 
may enter the dot at a much higher field than the other. 
It then moves towards its centre, delaying and possibly 
preventing the entry of a second main vortex. This ends 
up in a Landau state (e.g. for tp = +40°, FIG.EtO- No- 
tice that the vortex may continuously change its shape 
into a Bloch wall, provided that the dot is long and 
thick enough [l2j, |l9j. Thus, an asymmetric feature (like 
shape) is necessary to explain the experimental observa- 
tion of one- wall/ vortex state^3|- Notice also that more 
vortices may appear depending on the detailed shape of 
the dot (FIG. Eb,g). FIG. Ufc (upper part) shows that 
these simulations reproduce experimental H nuc (tp) con- 
vincingly. As a further step some simulations were per- 
formed in the finite-clement scheme to avoid numerical 
roughness on the surfaces. This yielded results quantita- 
tively similar to the case of vertical facets. 

Simulations of H ann (tp) (lower part of FIG.^) also re- 
produce experimental data. The jumps of H ann (tp) again 
result from the bifurcation before nucleation, implying 
different states before annihilation. The occurrence of 
one versus two H ann may be associated with the occur- 
rence of different flux-closure states at low field. A more 
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general experiment consists in proceeding to nucleation 
with the field decreased along a given angle ipi, followed 
by annihilation with the field increased along a different 
direction ip 2 . The plot H ann (if2), measured while keep- 
ing ifi fixed i.e. trying to prepare the system always in 
the same starting state, can be viewed as a signature of 
this state. FIG. 0Ji displays such plots for several val- 
ues of ipi, chosen in different branches of the experimen- 
tal H nuc (FIG. Et). The plots are not identical, which 
confirms that the remanent state depends on the angle of 
H cxt . The fact that different plots roughly consist of dif- 
ferent parts of a common set of two branches is also easily 
understood. For a given tp 2 a wall or a vortex of given 
circulation will always be pushed towards the same locus 
of the dot, be it alone at rcmancncc (vortex or Landau 
state) or having a companion (diamond state). In the 
latter case one point is found on each branch, whereas in 
the former case only one branch is revealed. 

Finally, the values of H sw i depend only weakly on 
the algorithm used (finite differences or finite elements). 
However, sometimes fine differences appear upon nucle- 
ation, such as the magnetization direction of vortices' 
cores, or the occurrence of more than two vortices (see 
FIG. c versus e; d versus f; g), although both codes 
were benchmarked successfully one against each other on 
time-resolved magnetization reversal issues [2^ ■ This un- 
derlines that describing the fine details of nucleation with 
simulations remains a challenge and that results should 
still be taken with care. 

To conclude, we reported the measurements, quanti- 
tative reproduction and understanding of angular nucle- 
ation and annihilation fields H mc (ip) and H^^ip) in a 
multi-domain magnetic particle. This is an example of a 
generalization to multi-domain states of the well-known 
Stoncr-Wohlfarth astroid. The main feature is the occur- 
rence of jumps in both plots, which result from the inter- 
play of a non-uniform magnetization state with the shape 
of the dot. Thus, such jumps should not be automatically 
ascribed to defects when observed in experiments. 

We thank Ph. David and V. Santonacci for technical 
support, and C. Meyer and W. Wulfhekcl for critically 
reading the manuscript. 
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FIG. 1: (a) AFM picture of the micro-SQUID. The dot 
strongly coupled to the micro-SQUID is indicated by a circle, 
(b) real-scale top view of the dot (c) Plot of H nuc (positive) 
and ffann (negative) for experiments (full symbols) and simu- 
lations (open symbols). Insets : experimental and simulated 
loops for two angles. 




FIG. 2: Schematic illustration of magnetization state bifurca- 
tion related to the orientation of the external field 7? ex t with 
respect to an edge. Arrows length sketch the magnitude of 
/f cxt . The external field H cxt is tilted (a) clockwise or (b) an- 
ticlockwise with respect to the normal to the edge. For a 
hexagonal dot, just after bifurcation the magnetization state 
may thus be in a leaf stateQ oriented either along the (c) long 
edge or (d) a diagonal. 
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FIG. 3: (a-d) Mid-height views of micromagnetic finite dif- 
ferences simulations of dots with vertical facets, for tp — 
40° (left) and <p = 50° (right), for (a-b) symmetric and (c- 
d) asymmetric in-plane shape. H ext is decreased from positive 
saturation towards the flux-closure state. Nucleation loci are 
indicated with full dots (e-f) Flux-closure state obtained on 
the same element with finite elements simulations (surface 
views) (g) Flux-closure states for ip = —40°, for finite differ- 
ences (left) and finite elements (right). The color codes the 
perpendicular component of magnetization (see scale). 
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FIG. 4: (a) Polar plots H iLnn (ip2) starting from different zero- 
field states, each state being prepared by initial saturation 
along a given direction ipi, shown in (b) with respect to the 
H nuc (tp) plot (c) Polar plot of // ann (yj) for radial field sweep- 
ing (same data as on FIG.|5J. 



